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Abstract. This paper presents the set of plane-parallel model atmosphere equations for a very hot neutron star (X-ray burst 
source). The model equations assume both hydrostatic and radiative equilibrium, and the equation of state of an ideal gas in local 
thermodynamic equilibrium (UTE). The equation of radiative transfer includes terms describing Compton scattering of photons 
on free electrons in fully relativistic thermal motion, for photon energies approaching c^. Model equations take into account 
many bound-free and free-free energy-dependent opacities of hydrogen, helium, and the iron ions, and also a dozen bound- 
' bound opacities for the highest ions of iron. We solve model equations by partial linearisation and the technique of variable 

I Eddington factors. Large grid of H-He-Fe model atmospheres of X-ray burst sources has been computed for 10' < T^ff < 3x10' 

, K, a wide range of surface gravity, and various iron abundances. We demonstrate that the spectra of X-ray bursters with 

^N) ' iron present in the accreting matter differ significantly from pure H-He spectra (published in an earlier paper), and also from 

blackbody spectra. Comptonized spectra with significant iron abundance are generally closer to blackbody spectra than spectra 
. of H-He atmospheres. The ratio of color to effective temperatures in our grid always remains in the range 1.2 < T^IT^ff < 1.85. 

' The present grid of model atmospheres and theoretical X-ray spectra will be used to determine the effective temperatures, radii 

fH ^ and MjR ratios of bursting neutron stars from observational data. 

^: 

^ ' Key words. Radiative transfer - Stars: atmospheres ~ Stars: neutron - X-rays: bursts ~ X-rays: spectroscopy 

1 . Introduction are recurrent events, but are not strictly periodic. The intervals 

between bursts are typically in the range of ~ 10"^ - 10^ s, but 

. ^ X-ray bursters were discovered by Grindlay et al. (1976), and sometimes a burster undergoes inactive phases lasting weeks 

^ , Belian et al. (1976). These X-ray sources are neutron stars in in- ^^^^ ^^^^^^ ^^^^ ^^^.^^ ^^.^^^^ j^^^ ^ Rappaport 1984; 

H ; teracting binaries. Very small mass functions and short orbital ^ewin, van Paradijs & Taam 1993). 

^ periods indicate that the secondary star has low mass (Stella rp, . . „ ^ u . ■ *u i a u ■ 

■ - - J V -pjjg origin or an X-ray burst is a thermonuclear flash in 

et al. 1987; Chakrabarty & Morgan 1998). Periodic pulsations ^ .i » j „ . * c t,,- •■ 

■' ° ^ ^ freshly accreted matter on the neutron-star surface. This idea 

have not been detected in the Hght curves of the great majority a u wr i e ^ /inn^\ a \/i u- e 

° J J was proposed by Woosley & Taam (1976) and Maraschi & 

of X-ray burst sources, and this has been used as an argument ^ ,. /irm\ t-.- j i i • n f.u 

° Cavaliere (1977). This model explains quite well many of the 

that the surface magnetic field of a neutron star is weak (Joss i u if » u . • i j- . 

. ° ^ global features of X-ray bursters, including energetics, tempo- 

& Li 1980) . We note that the spectrum of a type I X-ray burst j ^ tuu ■ f v u» 

' ^ -^^ ral structure and spectral behavior of an average X-ray burst 

becomes softer during the decay of the burst, and that such a imi imo a i- o i inoo\ 

° ^ source (Joss 1977, 1978; Ayasli & Joss 1982). 

softening is the signature of this type of events. . . • • ■ n j i ^ 

° ° -^^ Neutron star spectra were initially assumed to be perfect 

X-ray bursts are very energetic events. The estimated en- ^lackbodies. This assumption, however, is not valid for any re- 

ergy released during a single burst is typically ~ 10^^ ergs. The ^^.^^.^ ^j^^^^ atmosphere. For the case at hand, significant dis- 

light curve of a burst is chai-acterized by a rise time of ~ Is ^^^^-^^^ ^^^^ ^ ^-^^^^ blackbody specti'um result from radia- 

and a decay time of ~ 3-lOOs. X-ray bursts from a given source ^^^^^^^ neutron-star atmosphere, because the 

~ , „ . 7 . radiative opacity is dominated by scattering on free electrons 

Send ojfpnnt requests to: A. Majczyna r j jo 

e-mail: ma jczyna(acafflk. edu.pl (Madej 1974; Czerny & Sztajno 1983). In this case the outgo- 

' Pulsations wereobservedinSAXJ1808.4-3658,butotherfeatures ing X-ray spectra are shifted towards higher energies and the 

(e.g. kHz QPO) suggested that the magnetic field is relatively weak, peak X-ray flux is shifted in the same direction both for co- 

Recently at least 4 accreting millisecond pulsars have been detected, herent (Thomson) and incoherent (Compton) scattering. There 

and all are objects with weak magnetic fields. are two reasons for this effect. First, the free-free absorption 
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coefficient in the atmosphere of a hot neutron star varies with 
frequency as (for the opacity by heavier elements this re- 
lation is flat), and this causes a redistribution of radiation to- 
wards higher energies. Second, the presence of photon scatter- 
ing causes an additional shift of the outgoing flux in the same 
direction. The latter efi'ect results from the shift of the ther- 
malization depth to deeper (hotter) layers in the atmosphere. In 
other words, in the scattering atmosphere of a neutron star, pho- 
tons are created by thermal processes below the photosphere, 
where the local temperature T is much higher than the effective 
temperature Tgff, and where the peak of the emission spectrum 
is shifted to higher energies. 

An exact description of Compton scattering (Sampson 
1959; Pomraning 1973; Madej 1989, 1991a; Psaltis & Lamb 
1997; Psaltis 2001) is very important for the accurate deter- 
mination of basic neutron star parameters in X-ray bursters by 
use of fitting procedures. For this reason many authors have 
computed models of neutron star atmospheres with various de- 
grees of sophistication in the treatment of Compton scatter- 
ing. The influence of Compton scattering on the atmospheres 
of X-ray bursters has been discussed in several papers (see 
e.g. Ebisuzaki, Hanawa & Sugimoto 1984; and Foster, Ross 
& Fabian 1986). Consistent models of the atmosphere with 
Compton scattering have been presented by London, Taam & 
Howard (1986), Ebisuzaki (1987), and Ross & Fabian (1993). 
In these papers the authors utihsed the Kompaneets approxi- 
mation for the description of Compton scattering. 

We also note the works by Babul & Paczyhski (1987), 
Joss & Melia (1987), Zavlin & Shibanov (1991a,b) and Pavlov 
et al. (1991), who computed semianalytic model atmospheres 
and X-ray spectra of bursting neutron stars with luminosities 
approaching the Eddington hmit. Most of these papers pre- 
dicted the basic qualitative features of Comptonised model at- 
mospheres and spectra, such as the shift between color and ef- 
fective (or blackbody) temperatures. Another important feature 
found by these authors is the existence of a flux excess at low 
X-ray energies in the Comptonized emergent spectra of models 
with luminosities approaching the Eddington limit. 

Recently, papers with more general descriptions of 
Compton scattering were published by Madej ( 1 99 1 a,b), Madej 
& Rozahska (2000a,b), Joss & Madej (2001), and Madej, 
Joss & Rozahska (2004). The equation of transfer assumed in 
the latter papers allows for a large energy exchange beetwen 
photons and electrons, and therefore enables one to compute 
Compton scattering of photons with large initial energies, as 
well as the corresponding redistribution of photon energies dur- 
ing scattering. However, in these works the linear polarization 
of radiation was not considered, and the difl'erential Compton 
scattering coefficient was averaged over scattering angles. 

The present paper presents model atmospheres and spectra 
that will be used in the future for the quantitative interpretation 
of the observed X-ray spectra of bursting neutron stars. This 
paper is organised in the following manner: Section 2 presents 
the basic assumptions and equations of the model atmosphere 
with Compton scattering taken into account. Section 3 presents 
our computations of model atmospheres containing iron and 
the resultant theoretical X-ray spectra, and compares the re- 
sults with the H-He models of our previous paper (Madej, Joss 



& Rozahska 2004). Section 4 presents a summary and conclu- 
sions. 

2. Our models - assumptions and equations 

2.1. Assumptions 

We formulate the set of model neutron-star atmosphere equa- 
tions assuming: 

- plane-parallel geometry for the atmosphere, 

- local thermodynamical equilibrium and the equation of 
state of an ideal gas, 

- radiative and hydrostatic equihbrium, 

- zero magnetic field, 

- no neutron- star rotation. 

The equation of transfer includes free-free and bound-free 

absorption from all ions of hydrogen, helium and iron, as 
well as Compton scattering terms. The latter allow us to trace 
scattering off electrons with relativistic thermal motion, and 
with initial photon energies exceeding the electron rest mass. 
Moreover, the opacity of 10 spectral lines of highly ionised iron 
was included in our calculations. All these lines belong to the 
fundamental series of helium-like and hydrogen-like iron, and 
they fall in the energy range 6.7 - 9.0 keV, where most of the ra- 
diative energy of an X-ray burster is emitted. The total number 
of considered lines was chosen arbitrarily. 

We are aware that the assumed chemical composition prob- 
ably does not correspond to the atmosphere of a real X-ray 
burster. The composition of both the freshly accreted matter 
and the material processed during thermonuclear flashes in the 
stellar envelope will be rich of CNO and a variety of heavier 
elements. In our models, however, iron represents the average 
"metal", serving as the contributor of free electrons and the 
strong additional continuum free-free opacities. In this way we 
investigate the influence of both of these agents on the model 
atmospheres and the Comptonisation of the outgoing spectra, 
no matter what specific heavy elements may be the source of 
these efl'ects in a real X-ray burster. We are then able to com- 
pare the resultant spectra with the spectra of pure H-He atmo- 
spheres (Madej, Joss & Rozahska 2004). 

Ionization equilibrium in the model was determined in the 
following manner. For each discrete standard optical depth 
level T, and the corresponding values of temperature T and gas 
pressure Pg, we have guessed the value of electron pressure 
and then determined trial populations of all possible ionization 
states by solving the set of Saha equations for each element. 
Ionization populations computed in this way determine new 
values of electron density Ng, electron pressure Pe, and the re- 
sulting gas pressure P'J = P'JiPe, T) = {Nat + Ne) kT (ideal gas 
approximation). This procedure was iterated until Pg obtained 
in the above manner approached Pg with a relative accuracy 
better tiian 10"^ 

The energy-dependent absorption Ky in our model atmo- 
sphere code, ATM21, is the sum of many agents. First, we al- 
ways compute bound-free opacities from the 9 lowest levels of 
hydrogen, the 30 lowest levels of neutral heUum, and the 10 
lowest levels of singly ionized helium. Second, in this work 
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bound-free opacities of all existing iron ions from the ground 
level were computed following Vemer & Yakovlev (1995), who 
have published formulae fitting opacities of arbitrary elements 
and ions with atomic number Z < 30, stored in the Opacity 
Project database, see Seaton (1987). Third, free-free opacities 
of all existing ions were computed using the standard equations 
from Mihalas (1978). 

General relativistic efl'ects, such as the gravitational redshift 
of radiation and the bending of photon trajectories in strong 
gravitional field were not included in our equations. See Madej, 
Joss & Rozaiiska (2004) for a discussion of some of the effects 
of general relativity on the emergent spectra. 

The atmosphere of a neutron star with J,,// ^ 10^ K is 
geometrically very thin, usually of the order of 1-10 m (see 
Madej 1991b), except cases of the lowest values, g, for gravita- 
tional acceleration at the neutron star surface, when the model 
atmosphere approaches the critical gravity gcr- As g approaches 
gcr, the acceleration exerted on the atmosphere by the radi- 
ation field balances the downward gravitational acceleration, 
and the atmosphere expands to infinity. For higher values of g, 
the ratio of the height of the atmosphere to the neutron- star ra- 
dius is much smaller than unity. In this case the assumption of 
plane-parallel structure for the atmosphere is fully acceptable. 
The existence of radiative equilibrium implies that the radiation 
field is the only means of energy transport, and both convection 
and thermal conductivity can be neglected. Magnetic fields in 
X-ray bursts sources are relatively weak, and are in the range 
10^ - 10'' G (Joss & Li (1980); Lewin, van Paradijs & Taam 
1995; Miller, Lamb & Psaltis 1998). Consequently, the mag- 
netic field is unable to modify the continuum and line opacities. 

Our model equations describe the transfer of radiation sub- 
ject to the constraints of hydrostatic and radiative equilibrium. 
We assume the equation of state of an ideal gas in local thermo- 
dynamic equilibrium (LTE), so that all occupation numbers of 
bound and free states, opacities and emissivities are the same 
as their thermal equilibrium values at the local temperature T 
and density p throughout the entire atmosphere. However, in 
our model atmospheres the equation of transfer is dominated 
by the Compton scattering terms. This means explicitly that 
the equation of transfer distinctly remains in non-LTE. 

2.2. Model equations 

The equation of transfer with Compton scattering terms is as 
follows: 



/.^^=<(l-e-'"'/*")(Bv-/v) 

- IvS ^ r dv'criv^v',n-n')i\ 

X>' 47r Jo 



(1) 



where By is the Planck function, and ly and Jy denote spe- 
cific intensity and mean intensity of radiation, respectively. 
Parameter Ky denotes the true absorption coefiicient without 
correction for induced emission, whereas Ky = K'y{l - e'''^^'^^) 
will be used later as tiie same coefficient corrected for induced 



emission. Variable z denotes geometrical depth, and fi is the 
cosine of the zenithal angle (Madej 1991a). The variable 



cr(v -> v', n • n') = (Ty 0(v, v', n ■ n') 



(2) 



denotes the differential cross section of Compton scattering, 
where the redistribution function <p{v,V,n ■ n') is normalised 
to unity: 



(3) 



There exists the universal relation 

o-(v ^ v', « • «') v^e-''"/*^ = cr{V ->v,n' ■ n) v'^g-''"'/*^ , (4) 

which results from the detailed balancing of Compton scatter- 
ing in thermodynamic equilibrium (cf. Pomraning 1973). 

The dififerential cross section cr(v v',n ■ n') has been 
computed here following Guilbert (1981). Unfortunately, the 
set of the corresponding equations is very complex and cannot 
be presented explicitly here. 

We can integrate (piv, v', n ■ n') over the solid angles w', and 
obtain the zeroth moment of the redistribution function: 



(l)(y, V ) = (p (p(v, V , n - n ). 



(5) 



The equation of radiative transfer, Eq. 1, can be transformed 
to a useful form using the following approximations. First, we 
replace the redistribution function ((>(v,v',n ■ n') by its zeroth 
moment <l)(y, v'). Second, the specific intensity of radiation ly 
in both stimulated scattering terms was replaced by the mean 
intensity Jy. Then, we define two new Compton redistribution 
functions: 



<1>2(V,V') =1-1- 





/i(y-y') 







(6) 



3>(v,y'). (7) 



Third, we arbitrarily assume that the angular integral of cr(y' 
V, n ■ n') ly in Eq. 1 can be approximated by the product of two 
corresponding angular integrals of cr(y' v,n-n') and ly' . Eq. 
1 changes to 



dJy 

H— = Ky{By-Iy)-ayI, 

poz 



'y { ^i{v,V)dv' 
Jo 

r*00 

y I Jy,<^2{V,V)dV 

Jo 



(8) 



The equation of transfer can be written on the monochro- 
matic optical depth scale, dTy = —{Ky + ay)pdz ■ After simple 
algebra we obtain 



dl 

12—^ - ly-SyBy-il- Sy) Jy + 

OTy 



(9) 



{l-ey)Jy I ^i{v,V)dV -{l-Sy) I J y, ^2{v,V) dV , 
Jo Jo 

where By = Ky/iKy + (Ty) is the ratio of true absorption to the 
total opacity coefficients, and Compton scattering is described 
by the redistribution functions <l>i(y, y') and <l>2(y, v'). 
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Following the standard approach we obtain the zeroth and 
first moments of the equation of transfer in the following form: 



dTy 

dHy 

dTy 



= Hy (10) 
= ey{Jy -By) + {\- Sy) Jy I t/v' <^l{V, V) 

Jo 



I I dV 

Jo 



{I -By) dv'(i>2(v,v')Jy 



(11) 



Combining both the above equations we obtain the second or- 
der differential equation of transfer 



did . 
jyv^ l'7v^(/v-/v)) = Sy(Jy-B^) 



(12) 



(l-Sy)J, 

(l-fiv) 



I dv' <I>l(v,v' 
Jo 

Jr»co 
dv' (^*2(v,v')J 




where T]y - (Ky + o-y)/{K + o-)std, fv = Ky/Jy is the variable 
Eddington factor, and t denotes the standard optical depth. The 
latter is the monochromatic optical depth computed at the fixed 
(standard) wavelength. 

Terms with an asterisk are calculated at the unknown tem- 
perature at which the equation of radiative equilibrium is pre- 
cisely fulfiled. Since we seek this temperature, we make a lin- 
earization of the three following functions, expanding them in 
a Taylor series with respect to temperature 



B:(r) = B,(r) + Ar|^j , 

0*(y,v') = 0,(v,y') + Ar(-^ 
0;(y,v') = C)2(v,y') + Ar(^j , 



(13) 
(14) 
(15) 



where AT" - T* - T . The second and higher order derivatives 
were neglected in the above expansions. 

2.3. Radiative equilibrium 

The equation of radiative equiUbrium requires that 



(16) 



Differentiating the above equation with respect to t, and using 
Eq-Elyislds the useful equation of radiative equilibrium 



(17) 



0=1 dv rfySy {Jy - B'y) 
Jo 

+ I dv{\-Sy)Jy I dV^\{v,V) 
Jo Jo 

- I dv rjy(\ - Ey) I dv' <S>*2{y,v')Jy' 
Jo Jo 

The above equation is linearised according to the procedure 
described in Section 2.2. After that we obtain the temperature 



corrections in the form: 

dv Syrjy ( Jy " By) + L(t) 



AT = 



j^^dveyl]y(dBy/dT)-L'iT) 

where functions L(t) and L'(t) are defined as follows: 



L{t) 



= I dviJy-By) I dv'<^i(v,v' 
Jo Jo 



)x 



(18) 



(19) 



r]y(l - Sy) - l]y'(l - By') 



L'(t) 



Jo 



dv 



?7v(l - By) J, 

-T]y(l - Sy) 



I dv' 
Jo 

Jo '""'"^ 



(20) 



2.4. Hydrostatic equilibrium 

Our model assumes also the equation of hydrostatic equilib- 
rium, since we investigate static atmospheres. The condition of 
hydrostatic equilibrium can be written in the form: 



dP _ dPg dPr _ 
dz dz dz 



(21) 



Here the total scalar pressure P is the sum of gas and radiation 
pressures, P = Pj -H . If we write the above equation on the 
standard optical depth scale t, and use the expression 



dP 
It 



r 4n r~ 
c Jo 



dv. 



(22) 



then we get the final form of the hydrostatic equlibrium equa- 
tion: 



dPg ^ g An 
dT (k + o-),„d 



4jt n 
— riyHy 
c Jo 



dv . 



(23) 



2.5. The equation of transfer and boundary conditions 

The final equation of transfer, which is adopted in numerical 
calculations, is of the form: 

^ ylv-^ifvJv)] = ey (Jv - By) 



l-ey)Jy I Oi (y, v') c/y' - ( 1 - fv) I Jy ^liv, v') 
Jo Jo 

..^-(l-..)/.J^ ^dV.il-e.)l Jy^dV 

T/yfy (Jv - By) dv + L{t) 

l^7]y€y(dB/dT)dV-L'(T) 



(24) 



The equation can be solved when both the upper and the lower 
boundary conditions are specified. 

The upper boundary condition is of the standard form. 



— ifyJy) = hyJyiO) (25) 

OTy 
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where hy - Hy(0)/Jy(0) is the surface scalar factor which 
has to be iterated simultaneusly with the variable Eddington 
factors fy (Mihalas 1978). 

The lower boundary condition is in fact the thermalisation 
condition: 



(fvJv) 



1 dBy 

3 dTy 



(26) 



A detailed description and transformations of the inner bound- 
ary condition to the useful form are given in the Appendix. 

3. Results 

We calculated a grid of 106 models. Each of them was com- 
puted on a mesh of 144 standard optical depth points (12 
points per decade), distributed from Ts,d - 10"^ to lO'*, and on 
1139 frequency (energy) points, ranging from 400 keV to 12 
eV. Effective temperatures of our neutron star atmospheres are 
Teff-\, 1.5, 2, 2.5, 3x10^ K, and surface gravities § range from 
the critical gravity gcr to 10'^ (cgs units). Chemical composi- 
tion includes H and He in the solar abundance Nhc/Nh - 0.1 1, 
and iron either NfcINh = 3.7 x 10"^ (solar value), or 100 times 
higher. The constancy of bolometric flux across each model 
atmosphere was better than 0. 1 %, and only at a few points lo- 
cated deep below the photosphere was the flux error worse, up 
to about 2.5%. 

3.1. Temperature structure of models 

The distribution of temperature in a model atmosphere in radia- 
tive equilibrium is determined by several different parameters, 
such as the effective temperature, surface gravity, and the par- 
ticular chemical composition. Assumption of either Compton 
or Thomson scattering also is of crucial importance. In or- 
der to demonstrate particular features in T(t) stratification we 
have arbitrarily chosen models with the effective temperature 
Teff = 2 X 10^. In ah such models, the gradient of temperature 
is very steep in deep optical depths. In higher layers, e.g. for 
tross < 0.1, the run of T{T]ioss) exhibits different behaviour 

In the case of models with Thomson scattering, tempera- 
ture decreases with decreasing tross to some minimum value, 
and above this minimum T can rise only slightly and the model 
atmosphere is isothermal. A model atmosphere can be isother- 
mal starting from layers where it is optically thin, and photons 
do not efficiently interact with gas. In the case of a model at- 
mosphere with Compton scattering the situation is qualitatively 
different, since Compton scattering allows for energy and mo- 
mentum exchange between photons and free electrons. This 
process causes heating of the scattering layers of the atmo- 
sphere, because energy is transfered from radiation created in a 
hot layer of thermalisation to colder electron gas near the sur- 
face. Such an effect explains the existence of a temperature rise 
in the outermost layers, where Compton scattering dominates 
and sources of absorption vanish due to extremely low density 
and complete ionisation of atoms. 

On the contrary, a model atmosphere computed under an 
assumption of coherent Thomson scattering exhibits a mono- 



tonic run of temperature T, decreasing outwards. The isother- 
mal outermost region spreads above some optical depth lo- 
cated higher than in the Compton scattering model (see Fig. 
0. Moreover, the isothermal zone exhibits a much lower tem- 
perature. This is because Thomson scattering implies zero en- 
ergy exchange between photons and electrons, consequently, 
it cannot extract thermal energy from hard radiation created 
in hot thermalisation layers. In particular, in the model of 
Teff = 2x 10^ K, log g - 14.5, and the solar iron abundance, the 
atmosphere becomes isothermal at the Rosseland optical depth 
tross ~ 0.3 for Compton scattering and 0.03 for Thomson scat- 
tering atmospheres. Boundary temperatures in both cases are 
To « 2.21 x 10^ K and 1.45 x 10^ K, respectively. 




log ^BOSS 

Fig. 1. Temperature distribution in atmospheres of T^ff = 2 x 10^ K 
and logg = 14.5, and various chemical compositions. Both Compton 
and Thomson scattering models are shown here. The solid line repre- 
sents an atmosphere with solar iron abundance NfJNn = 3.7 x 10"' 
and the dashed line represents a model with the iron abundance 100 
times higher (Compton scattering models). Dotted line: model with 
solar iron abundance, long dashed line - 100 times higher iron abun- 
dance (Thomson scattering). 



In a model atmosphere with nonzero iron abundance, ther- 
mal absorption is higher than in a H-He atmosphere of the same 
Teff and logg. This is because iron is a strong opacity source. 
In particular, iron adds huge free-free absorption to the model 
atmosphere. Such an atmosphere is cooler in its uppermost lay- 
ers, see Fig. |2 This is because thermal nongrey absorption in 
LTE is a very efficient cooling mechanism due to the Kirchoff 
law, cf. also Mihalas (1978). Thermal emissivity and absorp- 
tion in LTE fulfil the equation jy = KyBy, where jy denotes 
emissivity, Ky is the thermal absorption coefficient, and By de- 
notes the Planck function. This relation implies that with in- 
creasing absorption, the rate of emission of radiation increases. 
If there exists a spectral window where the optical depth is low, 
then this layer easily loses thermal energy and its temperature 
T decreases. An atmosphere gets isothermal starting from the 
layer in which photons very weakly interact with matter. 
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Fig. 2. Temperature stratification in atmospheres of T^ff - 2 x 10' 
K and logg = 14.5, and different iron abundances. Solid line - solar 
abundance of iron, Nf^/N^ = 3.7 x 10"^, dotted line - iron abundance 
is 100 times higher; dashed line - H-He atmosphere. 

Surface gravity has an essential impact on the structure of 
an atmosphere (see Fig.|3}- In the case of high surface gravity, 
the density and thermal absorption coefficient Ky also are high. 
Consequently, photons are more efficiently absorbed and there 
exists a deeper minimum of temperature T than in the atmo- 
sphere of a neutron star with lower gravity. In a model with 
low surface gravity, the temperature T in the outermost lay- 
ers is higher than in a model with high gravity, since Compton 
scattering is the dominant opacity source and gas is effectively 
heated by multiple photon-electron energy exchanges. 
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Fig. 3. Temperature structure of atmospheres at the fixed Teff = 2 x 
10' K and solar iron abundance, and different surface gravities. Solid 
line - logg = 14.1, dashed line - 14.3, dotted line - 14.5, dot-dashed 
line - 14.9. 

The run of temperature in our models can be quite compli- 
cated in some cases, as is shown in Fig.0] The relation T(tross) 



displayed there is implicitly determined by a very complicated 
distribution of various iron ions and their continuum opacities, 
see also Fig.|5] 
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Fig. 4. Temperature structure of atmospheres with solar iron abun- 
dance and logg = 14.5, and of different r^yj; solid line - model of 
T,ff = 10' K, dotted line - 1.5 x 10' K. 



T 




Fig. 5. Fractional abundances for various iron ions in the atmosphere 
with Tfff = 10' K, logg = 14.5, and the solar iron abundance. 



3.2. Model spectra 

Theoretical X-ray spectra of hot neutron stars differ signifi- 
cantly from the blackbody spectra. In general, spectra of at- 
mospheres in which Thomson electron scattering opacity is 
present are harder than blackbody spectra (Madej 1974; van 
Paradijs 1982). This effect is also displayed in Compton scat- 
tering atmospheres, but then the frequency of the peak flux is 
lower. It is located between the maxima of a blackbody with the 
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effective temperature Tejf and the spectrum of an atmosphere 
with coherent Thomson scattering (Madej 1991a; Madej, Joss 
& Rozanska 2004). 

This effect is caused by two factors. First, adding of scatter- 
ing opacity to thermal absorption causes the layer of thermal- 
isation (where outgoing photons are created by thermal emis- 
sion) move to deeper and hotter layers as compared to a purely 
absorbing atmosphere. Consequently, the energy of the peak 
flux of the outgoing radiation increases. Second, noncoherent 
Compton scattering of hot photons in cooler surface layers re- 
moves (on average) part of the radiation energy. Due to this ef- 
fect, outgoing spectra are softer than those with Thomson scat- 
tering, and the previous increase of the peak flux energy must 
be reduced. 

Such a situation can be reversed by simultaneous redistri- 
bution of radiation flux caused by nongrey thermal absorption. 
In hydrogen-helium atmospheres of hot neutron stars both ele- 
ments are almost completely ionised, and then the only impor- 
tant thermal absorption is free-free absorption with Kv ~ v^^. 
Absorption strongly decreases with increasing energy, and it 
contributes to the additional increase of energy at the peak. 
However, the situation can be more complicated in atmo- 
spheres with nonzero iron abundance. 

As we can see in Fig. |6j energy of the peak flux of the 
outgoing continuum radiation of an atmosphere with Compton 
scattering is harder than in the case of the model atmosphere 
with Thomson scattering. In model atmospheres with iron, 
properties of outgoing radiation and its spectrum are deter- 
mined by the temperature structure in the atmosphere, which 
differs from T{t) in H-He atmospheres. In iron rich models, at- 
mospheres with Compton scattering are much hotter in opticaly 
thin layers than in Thomson scattering case (Fig.QJ. This jus- 
tifies the result that spectra with Compton scattering can some- 
times be harder than spectra with Thomson scattering (Fig.|6j. 

Comptonized spectra of hot neutron star atmospheres differ 
significantly from blackbody spectra of the effective tempera- 
ture Teff, both in the case of pure H-He and H-He-Fe chem- 
ical composition. The former model spectra were taken from 
Madej, Joss & Rozanska (2004). Differences between model 
atmospheres and blackbody spectra depend on the Tgjf, surface 
gravity log g and the chemical composition. 

All our spectra are harder than the blackbody of Te/f. 
Continuum outgoing spectra can be ordered in the following 
series, with increasing hardeness: the blackbody, spectrum of 
an atmosphere with iron abundance 100 times higher than the 
solar value, spectra of atmospheres with solar iron abundance 
{NfcINh = 3.7 X 10"^), and spectra of pure H-He atmospheres 
(see Figs.Q-|S|l. This is because in an atmosphere with low or 
zero iron abundance, the ratio of thermal absorption to the to- 
tal opacity coefficient is generally lower than in an atmosphere 
of higher iron abundance. In a series of atmospheres with the 
highest iron abundance NfcINh - 3.7 x lO"-', monochro- 
matic thermal absorption generally increases due to the pres- 
ence of strong bound-free continuum absorption of lithium-like 
to hydrogen-like ions of iron, and that absorption is strongest 
just beyond the corresponding spectral edge. Moreover, ionised 
iron causes a very strong free-free absorption, which does not 
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Fig. 6. Comparison of outgoing spectra computed either with 
Thomson scattering or Compton scattering, and the blackbody spec- 
trum. All model spectra displayed here are computed with T^ff = 
2 X 10' K and log g = 14.5. 

produce spectral features and therefore remains invisible in out- 
going spectra. 




2.5 3 3.5 4 4.5 



log ev 

Fig. 7. Outgoing flux spectra of the atmosphere of T^ff = 2 x 10^ K 
and logg = 14.5, and different iron abundances. Long dashed line - 
model with Nf^INh = 3.7 x 10"'', solid line - solar iron abundance 
NpelNn = 3.7 X 10"', dotted line ~ H-He atmosphere, dashed line - 
blackbody spectrum. 

Fig.|5]presents the sequence of outgoing theoretical spectra 
corresponding to the fixed T^jj = 2 x 10^ K and solar iron 
abundance NfcINh - 3.7 x 10"^, and various log g. We can see 
how theoretical X-ray spectra evolve towards higher energies, 
while the surface gravity log^ decreases towards the critical 
gravity gcr- 
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Surface gravity log^ significantly modifies the structure 
of an atmosphere because it directly determines stratification 
of density p through the equation of hydrostatic equilibrium. 
The spectrum of outgoing radiation emitted from the atmo- 
sphere of a neutron star of low surface gravity log^ is much 
harder than spectra corresponding to high log^. Atmospheres 
of low log^ are more extended and have lower densities than 
the atmosphere of a neutron star with higher gravity. Therefore, 
Compton scattering opacity gives a higher contribution to the 
total opacity in the former cases. The spectrum of an atmo- 
sphere with high surface gravity shows two prominent absorp- 
tion edges, caused by hydrogen-like and helium-like iron. We 
also see deep absorption lines belonging to the fundamental se- 
ries of both iron ions. 

We have included in our model atmosphere computations 
the series of lines 1 '5o - n^/^j of helium-Uke iron, and the se- 
ries 1^5 1/2 - of hydrogen-like iron. In both cases the 
principal quantum number of the upper level « = 2, . . . , 6, and 
higher Unes in both series were not included. Therefore we cal- 
culated profiles of a total of 10 lines. In this work we ignored 
the series of intercombination lines 1 '5o - rr'P^ of helium-like 
iron, which should be of noticeable strength in highly ionized 
iron. 

All the above lines cause a blanketing effect on model at- 
mospheres computed, and therefore they influence the radia- 
tive equilibrium and the run of temperature T{t) in all mod- 
els. Therefore, broadening of these lines is of particular impor- 
tance. We have considered iron lines broadened by natural (ra- 
diative), thermal (Doppler), and pressure (Stark) broadenings, 
the latter considered following Griem (1974). The final opac- 
ity profile of each of those lines was carefully computed as the 
convolution of all of the three mechanisms (Madej 1989). 




log ev 

Fig. 8. Flux spectra for T^ff = WK and logg = 12.9, and different 

chemical composition. Solid line - model with iron abundance 100 
times solar value, dashed line - the blackbody spectrum, dotted line - 
model spectrum of pure H-He atmosphere. 
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Fig. 9. Comparison of theoretical spectra of atmospheres with differ- 
ent logg and fixed T^ff = 2 x 10' K and solar iron abundance. Solid 
line - logg = 14.1 (it is very close to the Eddington limit), dotted line 
- 14.3, dashed line - 14.5, long dashed line - 14.7, and the dot-dashed 
line - 14.9. 

3.3. The ratio of color to effective temperatures 

Tables 1 and 2 present the ratios of color and effective tem- 
peratures for all the computed models. In the case of X-ray 
spectra computed with iron of solar abundance, the ratios of 
color and effective temperature do not differ strongly from 
those of pure H-He X-ray spectra (Madej, Joss & Rozanska 
2004). In both series of models (H-He atmospheres and those 
with solar iron abundance) we note the existence of local min- 
ima of relations TdTeff vs. logg for low effective tempera- 
tures (Jeff = 1 , 1 .5 X 10^ K). For higher effective temperatures 
{Tgff = 2 X 10^ or even higher) the ratio TdTeff increases with 
decreasing surface gravity. 

Table 1 shows that for low effective temperatures these 
ratios decrease with decreasing surface gravity from 1.35 (at 
logg = 15.0) to 1.28 (at logg = 14.3). The local minimum 
of TdTeff is rather flat, and its minimum value of 1.28 cor- 
responds to a wide range of surface gravities, ranging from 
logg = 14.4 to logg = 13.8. For lower logg the ratio 
TdTeff increases from 1.29 (for logg = 13.7) to 1.62 (for 
logg = 12.9), the latter value corresponding to a model which 
is very close to the critical surface gravity. The local minimum 
of TdTeff mentioned above separates regions on the logg axis 
where the outgoing spectrum (and the energy of peak flux) is 
determined mostly by the redistribution of radiation by highly 
nongray absorption and effects of Compton scattering, respec- 
tively. 

In spectra of neutron stars of high effective temperatures, 
the ratios TdTeff increase with decreasing surface gravity 
logg. This is because in very hot atmospheres the ratio of 
Compton scattering opacity to the total opacity coefficient at 
some fixed frequency increases and their spectra differ more 
and more from a blackbody. The highest value of the ratio 
of color and effective temperature in our iron-rich models ap- 



A. Majczyna et al.: Model atmospheres and X-ray spectra of iron-rich bursting neutron stars 



9 



proaches 1.77, and is obtained for the highest effective temper- 
ature of our grid, Teff = 3 X lO'' K, and log g = 14.8. We obtain 
quite a large value of TdT^fj - 2.01 also for log^ - 14.4, 
Tgff = 2.5 X 10^, and the solar abundance of iron. In the latter 
model 96% of the gravitational acceleration is compensated for 
by the acceleration exerted by the radiation field. 

For spectra of atmospheres with an iron abundance 100 
times higher than the solar value, the TdT^ff ratios are lower 
than those for atmospheres with lower iron abundance for the 
same Tgff and log g, see Table 2. This is because in atmospheres 
with overabundant iron the importance of thermal absorption 
is higher than in atmospheres with solar iron abundance. In 
models of low effective temperatures {Tgff = 1 x 10^ K and 
1.5 X 10^ K), the ratio TdTgff decreases with decreasing sur- 
face gravity log g. For Tgff = 2 x 10^ K and higher, these ratios 
increase with decreasing \ogg due to increasing significance of 
Compton scattering. 
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4. Summary and discussion 

We have computed and presented two grids of model atmo- 
spheres and theoretical spectra for hot neutron stars, which cor- 
respond to X-ray burst sources. The assumed chemical com- 
positions include hydrogen and helium in solar proportion, 
NhcINh = 0.11, and iron with two diff'erent abundances, 
NfcINh = 3.7 X 10"^ and 3.7 x 10"^ (the solar abundance and 
100 times solar). All spectra were calculated on a large grid of 
efl'ective temperatures T^ff and surface gravities log g. We cal- 
culated two grids of models bf with various iron abundances 
for 5 effective temperatures: 1, 1.5, 2, 2.5, and 3x10^ K, and 
with many surface gravities ranging from \ogg - 15 (in cgs 
units) down to the logarithm of critical gravity, log gcr, in steps 
of Alogg = 0.1. 

We assumed an LTE equation of state for an ideal gas, both 
the radiative and hydrostatic equilibrium, zero magnetic field. 
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Fig. 11. Run of ratios T^T^ff vs. log g in models with iron abundance 
100 times higher than the solar value. 

no neutron-star rotation, and plane -parallel geometry for the 
model atmospheres. The equation of radiative transfer includes 
free-free absorption from all ions, bound-free absorption from 
all ground levels of H, He, and Fe ions, and terms describing 
Compton scattering by free electrons. Bound-bound opacities 
(spectral lines) of helium-Uke and hydrogen-like iron are also 
included (10 lines of the fundamental series for both ions). The 
presence of scattering terms impUes that radiative transfer is a 
non-LTE process. 

All model atmospheres computed in this work exhibit a 
temperature inversion in the uppermost layers. This is caused 
by the heating effect of Compton scattering on the cooler elec- 
tron gas by hotter radiation from deep thermalisation layers. 
The same heating effect has been found in pure H-He model 
atmospheres (Madej, Joss & Rozanska 2004). 

Our calculations show that all model spectra of X-ray 
bursters diff^er from blackbody spectra at any given effective 
temperature Teff. The spectral differences depend on the pa- 
rameters of a particular model atmosphere, i.e., Tgff, log^, and 
chemical composition. The continuum spectra of iron-rich at- 
mospheres that were presented in this paper are always harder 
than a blackbody spectrum at a given T^ff. Similar results were 
also presented by Madej, Joss & Rozafiska (2004) for H-He 
atmospheres. Therefore, the color temperature Tc determined 
from the peak flux frequency is always higher than the effec- 
tive temperature Teff in neutron-star atmospheres dominated 
by Compton scattering. The ratio TdTgff does not exceed 1.85 
in our models. 

The theoretical spectra presented in this paper display two 
series of iron spectral lines, belonging to hehum-like and 
hydrogen-like iron (fundamental series). Line profiles are com- 
puted in LTE, which means that the Une source function is the 
Planck function. Profiles of these lines are computed as the nu- 
merical convolution of natural, thermal and pressure broaden- 
ing, the latter approximated on the basis of Griem (1976); see 
also Madej (1989). A few of the strongest absorption lines ex- 
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hibit reverse emission in the line cores, which was caused by 
the temperature inversion in the higher layers of the model at- 
mospheres. 

It is particularly interesting and new in our calculations that 
the peak fluxes of spectra with Compton scattering sometimes 
exhibit higher color temperatures Tc than Thomson scattering 
spectra computed for the same J,,// and log g. In the spectra of 
pure H-He atmospheres the situation is always different. In the 
latter case the peak fluxes of spectra with Compton scattering 
at a given Tefj are located at energies above the maximum for 
a blackbody spectrum but below the maximum for a Thomson 
scattering spectrum. 

This paper is an improvement and extention of the work 
by London, Taam & Howard (1986). These authors computed 
17 models of X-ray burst sources, but only 3 models assumed 
nonzero iron abundance (with a solar value of NpelNn - 10"^). 
Therefore, a comparison between their models and ours is 
very difficult. The TdT^ff ratios obtained by London, Taam 
& Howard (1986) are larger than our values. For example: for 
Teff = 1.17 X 10'^ K and log^ = 15.0 we obtained TdTeff = 
1.3, while London, Taam & Howard (1986) obtained 1.45; for 
Teff = 1.75 X 10^ K and \ogg = 15.0: 1.25 and 1.36, respec- 
tively; and for Teff = 1.25 X 10^ K and log^ = 14.0: 1.32 and 
1.46. 

Kuulkers et al. (2002) argued that if a burst aproaches the 
Eddington limit, then the deviation of its spectrum from black- 
body is the largest, TdTeff ~ 2, and that higher values of 
TdTeff do not occur. We confirm the claim with our precise 
numerical calculations. 

We note that our extensive grids of models are an excellent 
tool for fitting the observed spectra of X-ray bursts by use of, 
for example, the xspec software. Such spectral fitting is most re- 
liable when the atmosphere of the neutron star can be regarded 
as static, remaining in both hydrostatic and radiative equilib- 
rium. These conditions are best fulfiled in the cooling phase 
of an X-ray burst and during the quiescent intervals between 
bursts. 
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Table 1. T^Teff ratios in models with Nh/Nhc = 0.11, and Table 2. TdTeff ratios in models with Nn/NHe = 0.11, a 
Npe/Nn = 10-5. Npe/Nn = lO'l 
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Appendix A: The lower boundary condition 



At the lower boundary Tsui - Tmax we assume that the radiation 
field is thermalised at each frequency of the grid, Jv(Tmax) - 

The second moment of the radiation field equals 



Neglecting terms with (AT)^ we obtain: 



Ky + CTv dT 



X 



(A.ll) 



1 +AT 



dT 



M2 

dT^ Ml 



Ky^fyJy^fyBl, (A.l) 

1 .1 T-jj- r T 1 /T 1 1 J Temperature corrections AT are necessary to complete the 

where the Eddmgton factors fy - 1/3. The asterisk attached , ^ . r„, . , , • , ^ 

^, „i 1 ^ ■ J- 1 . ■ 1 ■ 1 ^ 1^:1 1 above equation. These corrections can be obtained from the 

to the Planck function indicates that Its value strictly fulfils the . ,. ... . 

... ^ ,. . .... . equation of radiative equilibrium: 

condition of radiative equilibrium. 

Differentiating of Eq. A.l yields: 

1 dBl dT 



dTy 3 dTy 3riy(K + o-)std P dT dz 



(A.2) 



where dXy - -T]y(K + o-)stdpdz, the ratio rjy - (Ky + o-y)/(K + 
cr)std , and the dimensionless absorption Sy = Ky/(Ky + cTy). 

The temperature gradient dT/dz, which is required in the 
above equation, can be determined from the bolometric flux: 



H = 



1 



dv 
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dBl dT 



(Ky + CTy) p dT dz 



Bolometric flux is given by: H - o-rT*, J4jt , where ctr 



5.66961 X 10 erg cm ^ s ' K Temperature gradient is 



dT 

— = 3//X 

dz 

therefore, 
_d_ 
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(fvJy) 
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dBl 



T]y(K + cr)„dp dT 

1 dBl 



(A.5) 
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(Ky + (Ty) P 8T 

We approximate dBl/dT in the above equation by a linear 
Taylor series 



dBl ^ dBy ^ d^By 
dT dT dT^ 



(A.6) 



and obtain the lower boundary condition in the following ex- 
panded form: 
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(A.l) 



f 

Jo 



dv ■ 



1 



;o (Ky + a-y)p\dT 
Let us introduce new variables: 

Ml = dv , 

Jo oT Ky + a-y 



dBy^^^d-B. 



dT^ 



M2 = 



dv 



d^By 



1 



dT^ Ky + (Ty 

The lower boundary condition assumes the following form 



(A.9) 



dT, 



(fyJv) 



HM 



dT dT^ 



(A. 10) 



/-•CO 

dv r]ySy(Jy - B*y) 

Jo 

r^oo /-•CO 
+ dv (l-Ey)Jy dv'^\(vy) 

Jo Jo 

/-•DO /-•OO 

- I dvT]y(l-ey) I dv (^l(v,v')Jy, . 
Jo Jo 



(A. 12) 



EQ. IA.12l should be linearized with respect to temperature T, in 
(A 3) ^'^ replace variables with an asterisk: 

r(Ty)^T(T,) + AT(Ty), 

ldBy\ 



Bl(T) = By(T) + AT 



(A.13) 
(A. 14) 



(A.4) (l)?(v,v') = <l>i(v,v') + Ar 



'dT 



cD*(v,y') = 02(y,v') + ArJ^j 



This step yields temperature corrections in the form: 

dv SyJ]y (Jy " By) + L(t) 



AT = 



dv eyr^y(dBy|dT)-L'(T) ' 
where functions L(t) and L'(t) are defined as: 

/-•OO /-•OO 

I dv(Jy-By) I dv'^i(v,v) 
Jo Jo 



L(t) 



')x 



(A.15) 
(A. 16) 

(A.17) 
(A. 18) 



r]y(l - Sy) -T]y'(l - Sy') { — 
V 



L'(T) 



-r 

Jo 



dv 



r]y(l - Sy)J, 

-?7v(l - Sy) 



Jo dT 



(A. 19) 



(A. 8) Combining equations A.ll and A.17 we obtain the final form 
of the lower boundary condition: 



^ T, -^"^r' dBy 
■7- (fvJv) = — ^ X 

dTy Ky + (Ty til 



(A.20) 



X O -I- 



Jy dVSy7]y( Jy- By) + L(t) 
1^°° dVSyl]y(dBy/dT)-L'(T)' 



X (Ml + M2AT) 



1 



dBy 

Ky + (Ty dT 



1 +AT 



dB, 



' d^By M2 



dT ] dT^ M, 



AT 



dBy 
dT 



d^By _ Ml 



